Massive compact dwarf stars and C-field 
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Abstract 

The effect of C-field in high density matter has been studied. We find that the 
negative energy and negative pressure of the C-field helps in formation of massive 
compact stable neutron stars of mass ~ 0.5 solar mass which is in the range of 
0.01 to 1.0 solar mass of recently observed dwarf stars. 
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There has recently been reported observations by gravitational microlensing of dark dwarf stars 
in the mass range of 0.01-1. OM (Alcock et al. 1993; Aubourg et al. 1993). Simple extrapolations 
of this observations have led people to speculate (Boughn & Uson 1995) that such bodies could 
be numerous and make up the "dark matter" inferred from the rates of rotation of galaxies 
and superclusters of galaxies. In any case, this observation has opened up the question of what 
these massive compact objects themselves are. Cottingham, Kalafatis & Vinh Mau (1994), have 
proposed a model where these objects are identified as quark stars formed after quark-hadron 
transition. 

We, in this work, suggest a different understanding of these objects. The stability of these 
objects at unusual mass values is indicative of a simultaneous increase in binding and reduction 



in the internal pressure compared to the normal stars leading to the balance between the two at 
lower mass values. Driven by such an argument and since the creation field (C-field) of Pryce, 
Hoyel & Narlikar (1962) and Narlikar (1973) used in the context of steady-state theory has 
exactly these characteristics of negative pressure and negative energy, we have tried to include 
its contribution into the description of high density matter interms of SU(2) chiral sigma model 
(Sahu, Basu & Datta 1993). 

Eventhough, the hot big-bang model of the creation and evolution of the Universe has gained 
acceptance over the steady-state cosmology, it does have problems associated with linearity of 
Hubble flow and determination of the age of the Universe from the Hubble Space Telescope data 
(Narlikar 1993). It is, therefore, in the fitness of things to explore the possibility of an explanation 
for the problem at hand within the context of an alternative model as has been attempted in 
other contexts (Weinberg 1972; Arp et al. 1990). 

The massless C-field was originally introduced by Pryce and later extensively used by Hoyle & 
Narlikar (1962) and Narlikar (1973) to provide a field-theoretic understanding of the continuous 
creation of matter in the steady-state evolution of the Universe preserving energy and momentum 
conservation. As a first approximation, we take the C-fields to be noninteracting. Its energy- 
momentum tensor is given by 



where Cj is simply |£ and f (> 0) is a coupling constant. In fact, the energy density and pressure 
contributed by C-field are given by 



T ik = _f^ C i C k _ _ g ikC l Ci\ 
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where C = ^E-. Again, if all particles, with whose creation the C-field are simultaneously created, 



have same mass m, then C = m. On quantisation of C-field, it has been argued that the C-field 
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does not lead to the usual problem of cascading associated with negative energies once its effect 
of expanding the space-time structure due to the feed back of energy-momentum tensor of C-field 
on space-time geometry is taken into account. However, as is borne out by our calculation, the 
negative energy of C-field makes matter to condense to more compact size of mass of the order 
reported above. 

We like here to state that the equation of state (EOS)of matter above the nuclear matter 
density plays a crucial role to determine the equilibrium structure of compact objects such 
as neutron stars. Upto about nuclear density the equation of state is reasonably well known, 
but the central density of compact objects can be almost an order of magnitude higher. In 
this regime, the physics is unclear. For that reason, attempt has been made to calculate EOS 
using two major techniques: (i) non-relativistic and (ii) relativistic field theories. The shape 
of baryonic potential is not known at very small interparticle separations. Also, it is not clear 
that the potential description will continue to remain valid at such short ranges. Therefore, the 
non-relativistic approach may not be adequate. In the recent times, the relativistic approach 
has drawn considerable attention. In the relativistic approach, one usually starts from a local, 
renormalizable field theory with baryon and explicit meson degrees of freedom. The theory is 
chosen to be renormalizable in order to fix the coupling constants and the mass parameters by 
empirical properties of nuclear matter at saturation. As a starting point, one chooses the mean 
field approximation which should be reasonably good at very high densities. This approach is 
currently used as a reasonable way of parameterizing the EOS. However, in recent years, the 
importance of the three-body forces in the EOS at high densities has been emphasized by several 
authors (Jackson, Rho & Krotscheck 1985; Ainsworth et al. 1987). This gives theoretical impetus 
to study the chiral sigma model, because the non-linear terms in the chiral sigma Lagrangian 
can give rise to the three-body forces. 
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We also take the approach that the isoscalar vector field be generated dynamically. Inclusion 
of such field is necessary to ensure the saturation property of nuclear matter. The effective 
nucleon mass then acquires a density dependence on both the scalar and the vector fields, and 
must be obtained self-consistently. We do this using mean-field theory wherein all the meson 
fields are replaced by their uniform expectation values. 

The Lagrangian for an SU (2) chiral sigma model that includes an isoscalar vector field {uj^) 
is (h — 1 = c) (Sahu, Basu & Datta 1993) and a non-interacting C-field is 

C = l -{d^.d^ + dpo&o) - ±{t.!t + a 2 - x 2 f 

-\f^f^ v + \gJ{° 2 + 1? 2 )^" 

—G^CT + \my^ ~ \gMlt r~*-f)1> ~ td.Cd^C, (3) 

where F^ v = d^u u — d u u^, G^ v = d^p u — d u p fM , ip is the nucleon isospin doublet, if is the 
pseudoscalar pion field and a is the scalar field. The expectation value < ■07 o '0 > is identifiable 
as the nucleon number density, which we denote by rig- 

The interactions of the scalar and the pseudoscalar mesons with the vector boson generates a 
mass for the latter spontaneously by the Higgs mechanism. The masses for the nucleon, the scalar 
meson and the vector meson are respectively given by m = g a x \ m a = \/2\x ; m u = g^x Q , 
where x Q is the vacuum expectation value of the sigma field. C-field being non-interacting remains 
massless. 

The equation of motion for the mean vector field specifies u Q 

uj =^- 2 , x = (<^ + ^ 2 >) 1/2 - (4) 
The equation of motion for a written for convenience in terms of y = x/x Q is of the form 
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where m* = ym is the effective mass of the nucleoli and c a = g^/m 2 ; c w = g^/rri^. 

At high densities typical of interior of neutron stars, the composition of matter is asymmetric 
nuclear matter with an admixture of electrons. The concentrations of protons and electrons can 
be determined using conditions of beta equilibrium and electrical charge neutrality. We include 
the interaction due to isospin triplet p- meson in Lagrangian for purpose of describing neutron-rich 
matter. The equation of motion for , in the mean field approximation , where ~7f is replaced 
by its uniform value pi (here superscript 3 stands for the third component in isospin space), 
gives pi = ( n p — n n)- The symmetric energy coefficient that follows from the semi-empirical 
nuclear mass formula is a sym = + 6(fc 2^ 2)l/2 , where c p = g 2 p /m 2 p and k F = (6tt 2 pb/7) 1/3 
(Pb — n P + n n) ■ We take the values of c a , c u and c p by fits of saturation density (0.153 fm~ 3 ), the 
binding energy (-16.3 MeV) and symmetric energy (32 MeV) (Moller et al., 1988) in the absence 
of C-field as a first approximation. These give c a = 6.20 fm 2 ; c w = 2.94 fm 2 ; c p = 4.6617 fm 2 . 

The diagonal components of the conserved total stress tensor corresponding to the Lagrangian 
together with the equation of motion for the fermion field (and a mean field approximation for 
the meson fields) provide the following identification for the total energy density (e) and pressure 
(P) for neutron star system : 



m 



2 (i-y 2 ) 2 i 2 cUK 3 + fcn 3 ) 2 7 2 c,fe 3 -fc n 3 ) 2 

8c CT 72n 2 y 2 72n A y 2 





P = 



2 (l-y 2 ) 2 ! 2 Uk p 3 + k n 3 ) 2 7 2 c.(fe P 3 - K 3 ) 2 
8c a 72n 2 y 2 72vrV 





A specification of the coupling constants c a , c u , c p 



and / now specifies the EOS. 



As of the coupling constant /, since it has the dimension of (mass) 2 , we have parameterized it 
in two ways like, / = fm 2 and / = /— , where m is typically the nucleon mass, tib is the baryon 
density and / is dimensionless. The second parameterization employing a linear proportionality 
between / and n# is more reasonable and physical as the coupling of C-field is likely to grow 
in strength with increase in % since creation of baryons occurs in association with creation of 
C-field. 

The structure of a neutron star is characterized by its gravitational mass (M) and radius (R). 
These gravitational mass and radius for non-rotating neutron star are obtained by integrating 
the structure equations, which describe the hydrostatic equilibrium of degenerate stars : (Misner, 
Thorne & Wheeler 1970) 

dp G(p + p/c 2 )(m + 47rr 3 p/c 2 ) dm , 9 

— - = — — — — — — — -; = Anr p, (8) 

dr r 2 (l — 2Gm/rc 2 ) dr 

where p and p(e/c 2 ) are the pressure and total mass energy density. For a given EOS, p(p), 
and a given central density p(r = 0) = p c , eqs.@ are integrated numerically with the boundary 
condition mir = 0) = to give R and M. The radius R is defined by the point where P ~ 

0, or, equivalently, p = p s , where p s is the density expected at the neutron star surface. The 
maximum total gravitational mass for stable configuration is then given by: M = m(R). The 
results obtained after the inclusion of contribution of the C-field are given in table I and figure 

1. In table I, we have presented the maximum mass (M) and the corresponding radius (R) for 
a more typical stable star as a function of central density (p c ). It is observed from the figure 
that with increase of /, the EOS becomes softer leading to reduction of stable neutron star 
gravitational mass. 

A comparison of the values of the radii (R) and gravitational mass (M) for various values of 
dimensionless coupling parameter / with those in the absence of C-field as obtained in an earlier 
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work (Sahu, Basu & Datta 1993) (case III of Table I) reveals that the size and mass of stable 
neutron star structures have been significantly reduced. This can be understood on the ground 
that the negative energy supplements gravitational binding and the negative pressure reduces 
opposition to gravitational agregation of matter. Thus, C-field doubly facilitates formation of 
compact stable neutron star structure of smaller mass and dimension. In fact, for / values near 
0.05 and 1.5 in the two schemes of parameterization of dimensional coupling constant /, the stable 
neutron star mass is about O.5M which is well inside the mass range of dwarf stars recently 
observed between 0.01 to 1.0M Q considered to be candidates for the dark matter ( Boughn & 
Uson 1995; Cottingham, Kalafatis & Vinh Mau 1994). Further, the value of / being of 0(1) 
for the mass of compact object to lie in the right range implies that the interaction involved is 
strong in character as is to be expected between the nucleon and scalar field. 

Apart from the above agreement, we believe that the study of the effect and role of the 
C-field in various astrophysical problems needs to be taken up in its own merit. This provides 
scope to investigate the existence of explanations alternative to those found within the ambit 
of the generally accepted big bang model of the Universe. The present work completes a small 
programme in that direction. 
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TABLES 

TABLE I. The maximum mass (M), the corresponding radius (R) and central density (p c ) of more 
typical compact dwarf stars for various values of dimensionless coupling parameter / for three cases; 
Case I: / = fm 2 , Case II: / = /ng/m and Case III: / = 0. 



Pc 


R 


M/M Q 


/ 


Cases 


(g cm -3 ) 


(km) 








2.0xl0 15 


10.96 


2.16 


0.001 




LOxlO 15 


7.61 


1.51 


0.005 


I 


7.0xl0 15 


5.94 


1.19 


0.01 




2.5xl0 16 


3.12 


0.64 


0.05 




3.0xl0 15 


9.10 


1.678 


0.5 




7.5xl0 15 


5.66 


1.01 


1.0 


II 


2.0xl0 16 


3.61 


0.67 


1.5 




1.5xl0 15 


13.65 


2.59 


0.00 


III 
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FIGURES 

FIG. 1. Pressure and energy density curves for three cases; Case I: / = fm 2 with / = 0.001 (curve 
b) and / = 0.05 (curve c); Case II: / = /ng/m with / = 0.5 (curve d) and / = 1.5 (curve e); Case III: 
/ = (curve a). 
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